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ABSTRACT 

Initial research on the external flow, low mass liquid sheet radiator (LSR), has been 
concentrated on understanding its fluid mechanics. The surface tension forces acting at the 
edges of the sheet produce a triangular planform for the radiating surface of width, W, and 
length, L. It has been experimentally verified that L/W agrees with the theoretical result, 
L/W = [We/8] 1/2 , where We is the Weber number. Instability can cause holes to form in 
regions of large curvature such as where the edge cylinders join the sheet of thickness, r. The 
W/t limit that will cause hole formation with subsequent destruction of the sheet has yet to 
be reached experimentally. 

Although experimental measurements of sheet emissivity have not yet been performed 
because of limited program scope, calculations of the emissivity and sheet lifetime as 
determined by evaporation losses were made for two silicon based oils; D.C. 705 and Mez- 
Emissivities greater than 0.75 are calculated for t > 200 pm for both oils. Lifetimes for Me^ 
are much longer than lifetimes for 705. Therefore, Mej is the more attractive working fluid 
for higher temperatures (T > 400 K). 


CURRENT STATUS OF LIQUID SHEET RADIATOR RESEARCH 
Introduction 

The liquid sheet radiator (LSR) is an external flow radiator similar in several ways to other 
external flow radiators such as the liquid droplet [1-3] (LDR) and liquid belt radiators [4,5]. 
All of these radiator concepts potentially have lower mass than solid wall radiators such as 
pumped loop and heat pipe radiators. They are also nearly immune to micrometeoroid 
penetration. However, the LSR has the added advantage of simplicity. As a result of surface 
tension forces acting at the edges of the sheet, the sheet flow has a triangular shape (Fig.l) 
which is desirable as it allows collection of the flow at a single point. To obtain a similar 
triangular area for the LDR requires very accurate aiming of many (10 5 to 10 6 ) droplet streams. 
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(a) Top view. 



(b) Cross section view. 


Figure 1 . — Schematic of thin liquid sheet flow. 


Simplicity of the LSR should also result in lower cost and mass for the LSR compared to the 
LDR and liquid belt radiators (LBR). 

A major problem for all external flow radiators is the requirement that the working fluid have 
a very low (<10' 8 Torr) vapor pressure to minimize evaporative losses. As a result, working 
fluids are limited to certain silicon based oils (such as used in diffusion pumps) for low 
temperatures (300 to 400 K) and liquid metals for higher temperatures [I]. 

Figure 1 is a schematic drawing of a liquid sheet flow. Surface tension forces at the two 
edges of the sheet push the edges toward the z axis. As a result, as the flow moves in the 
z direction the edge cross-sectioned area, grows. In order to satisfy mass continuity the 
edges approach each other and finally meet at the point z = L. Figure 2 is a photograph of 
a small scale sheet flow of Dow-Corning 704 oil in vacuum. The growth of the edge cylinders 
can be seen in this photograph. If a gas surrounded the sheet flow the associated drag force 
would tend to break up the sheet flow. Whether or not sheet breakup occurs depends on the 
velocity and ratio of the sheet liquid density to the surrounding gas density. 

Currently the LSR is only at the research stage of development. No system hardware, other 
than the narrow slits used to form the sheets, has been developed. However, a schematic 
drawing of a possible LSR space radiator is shown in Fig. 3. This shows four sheet flows 
connected in parallel. The number of sheets required for an actual application will be 
determined by the heat load. The system consists of three main components; a slit ejector, a 
collector, and a pump to produce the necessary flow. In Fig. 3 a heat exchanger is shown 
between the thermal load and the sheet flows. This heat exchanger is not necessary if the 
thermal load is compatible with the sheet liquid. The temperature drop that occurs between 
the ejector and collector is very small so that the LSR behaves like a constant tempertature 
radiator, similar to a heat pipe radiator. 

Initial research on the LSR has been directed at understanding the fluid mechanics of sheet 
flows [6,7]. There has also been a systems study [8] of the LSR in low power (2 kW) Closed 
Brayton (CBC) and Free Piston Stirling (FPS) cycles. In that study, with the LSR mass 
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Thermal load 



Figure 2. — Sheet flow of Dow-Corning 704 oil 
in vacuum. Slit thickness = 100 jam, slit 
width = 3.4 cm, L/W = 9.6. 



Figure 3. — Schematic of parallel connected liquid 
sheet radiator system. 


estimated to be 1,5 kg/m 2 the LSR was shown to be quite compatible with FPS based space 
power systems, but not well matched with heat rejection temperature profiles of CBC power 
systems. However, general conclusions from the study confirmed the applicablility of LSR to 
space power systems with near constant heat rejection temperatures, including Stirling, 
Rankine, Thermoelectric, and Thermionic power systems. 

This paper presents the latest results for the scaling (L/W = sheet length/sheet width) and the 
geometry of the edge cylinders. Also, since emissivity and evaporative mass loss are the two 
most critical properties of a sheet radiator, results of an emissivity calculation and sheet 
lifetime results based on experimental evaporative loss rates for two possible oils that can be 
used in the LSR are presented. Finally, mention of LSR research currently in progress and 
important research yet to be done will be made. 


Scaling Results 

The surface tension force, dF, at the two edges of the sheet push the edges toward the z 
axis (Fig. 1). This force is proportional to the surface tension, a, and inversely proportional 
to the radius of curvature Re- As the flow moves in the z direction the edge cross-sectional 
area, A,., grows. In order to satisfy mass continuity the edges approach each other and finally 
meet at the point z = L. 
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If the surface tension force per unit length, is equated to the rate of momentum change 

dz 

per unit length, pm/, where p is the density, t is the sheet thickness, and u e is the edge 
velocity, then relations for u e and LAV can be obtained [7]. 



For the case where there is no force (gravity for example) acting in the z direction the sheet 
thickness r, is a constant. As a result the edge velocity, u e , will be constant [7] and; 



Where the Weber number, We, is defined as: 


We = 


pwj 

a 


( 2 ) 


(3) 


where w 0 is the velocity in the z direction. As can be seen from Eqs. (2) and (3), LAV is 
a linear function of velocity, w D . 

A large amount of L/W data for D.C. 705 oil has been obtained in a vacuum facility [7] 
which can accommodate sheet flows with W ^ 25 cm and L < 3.5 m. These results are shown 
in Fig. 4. Note the good agreement between the theoretical result (Eq. (2)) and the 
experimental data. For a practical radiator L/W ® 10 so that We « 1000. The temperature 
dependence of p and cr for Dow-Coming 705 is given by [9]: 

p = 1569 - 1.625 T kg/m 3 T = K < 4 ) 


a 


86.14 - 



10' 3 n/m 


T = K 


(5) 


All the data in Fig. 4 was taken with T = 371 K. For D.C. 705 oil 0.025 < <r < 0.040 n/m 
for 400 > T > 300 K, which is the temperature range where the silicone oils have sufficiently 
low vapor pressure. For higher temperature operation liquid metals [1] are the only fluids with 
low enough vapor pressure to have minimal evaporative losses for space applications. The 
surface tension [1], a, of liquid metals is at least an order of magnitude larger (0.3 to 0.9 n/m) 
than that for silicone oils. Whereas, the densities of the liquid metals of interest range from 
one-half (lithium) the density of D.C. 705 to 12 times D.C. 705 density for mercuiy. 
Therefore, for a given sheet thickness, t, the low density metals such as lithium and aluminum 
require a significantly higher velocity than a silicone oil to obtain the same Weber number. 
For a silicone oil w c » 10 m/sec yields We * 1000 for t = 200 pm. However, for a low 
mass liquid metal such as lithium w 0 « 60 m/sec for We = 1000 and t = 200 pm. The 
other liquid metals will require velocities of 10 to 60 m/sec to obtain We = 1000 for 
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Figure 4. — Comparison of theoretical and experimental sheet scaling. 


t = 200 |irn. The higher velocities required by the liquid metals coupled with their higher 
densities translates to significantly larger pump power than that required for the silicone oils. 

Sheet Stability and Edge Cylinder Geometry 

Another critical fluid mechanical problem is sheet stability. Linear stability theory for an 
infinite sheet [11] predicts no instability (i.e., no solutions that grow in time). However, if 
the sheet is not planar everywhere, such as occurs at the sheet edge cylinders, then linear 
theory does predict growing solutions in the nonplanar regions. These results were obtained 
in a study that has not yet been completed and hence will be reported elsewhere. However, 
some highlights and general observations are discussed below. 

Holes have been observed to form where the edge cylinder joins the sheet but are quickly 
swept downstream before they have a chance to grow significantly. However, as x decreases 
the hole growth rate will increase. It is expected that there will be a limit on W/x to insure 
that a hole will not form and then grow fast enough to eliminate the sheet. Experimentally this 
W/x limit has not been reached. However, an experiment is being set up to determine it. 

Since instability is expected to occur in the edge cylinder region, a study of the edge cylinder 
geometry was made. In Fig. 5 a representation of the edge cylinder geometry at four z (flow 
direction) locations is shown. Experimentally it has been found that the edge cylinder 
geometry oscillates in the z direction between nearly a circle to a highly elliptical shape then 
a "peanut" like shape and finally back to nearly a circle as shown in Fig. 5. 
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(a) Initial, nearly circular cross-sectional shape. 


(b) Flattened out cross-sectional shape. 



Dimensionless edge cylinder area, A c = Wz/tL 


Figure 6. — Comparison of experimental and theoretical edge cross section dimensions. 


6 





The theoretical analysis predicts a nearly circular geometry for the edge cylinder. No highly 
elliptical or peanut like solutions were found. In Fig. 6 a comparison between the experimental 
and theoretical edge cylinder geometry results is made. Experimental results for three sheet 
thicknesses, x = 50, 80, and 100 pm are shown. The dimensionless width of the edge 
cylinder, r = 2r/x, and dimensionless maximum height of the edge cylinder, Hmax = 2s/x, are 
plotted as a function of the dimensionless edge cylinder cross sectional area, 
A c = 2A c /f 2 = Wz/tL (from mass continuity and the fact that the sheet is triangular it can be 
shown that the edge cylinder area is a linear function of z)._The results for the x = 100 pm 
sheet agree well with the model shown in Fig. 5. For small A c the experimentalr and rj MAX 
values are close to the nearly circular theoretical results. (For a perfect circle r = 2 Hmax ) 
However, as increases, r increases while r| MAX decreases below the theoretical result, 
indicating the highly elliptical shape is being formed. Then at A, » 75, r starts to decrease 
toward the theoretical result as r| MAX increases toward its theoretical result, indicating that a 
nearly circular shape is again being formed. Then, at \ a 100 the cycle begins again. Only 
from photographic results is it known that the peanut like shape occurs. For the x =_50 and 
80 pm cases both r and r| MAX remain above the theoretical result for all A,.. For 
x = 80 pm a single oscillation from nearly circular geometry to elliptical shape and then back 
to a circular shape can be seen. Whereas for x = 50 pm two^such oscillations are seen. 

All the data in Fig. 6 indicate an elliptical geometry with r > 2 r) MAX . The more flattened 
the cross section (r » 2 r| MAX ), the smaller the curvature (l/RJ will be at the point where the 
edge joins the sheet. The instability growth rate is proportional to the curvature. Therefore, 
the flattened edge cross section regions are less likely than the more circular regions to result 
in hole formation. 


Properties of Possible Silicone Oils for the LSR 

Low vapor pressure is a necessity for the working fluid of an external flow radiator. 
However, for an efficient radiator the working fluid emissivity must also be large. Fortuitously, 
for the temperature range of interest (300 to 400 K) we find that the silicone oils, which are 
transparent to visible radiation, do have sufficiently high emissivity for space radiator 
applications. This emissivity results from a high absorption coefficient in the wavelength 
region (8 < X < 15 pm) where most of the radiation occurs for temperatures of 300 to 400 K. 
In other words, the silicone oils behave like a black body in the 8 < X < 15 pm region, as 
shown in Table 1. 

Previously [7] the emissivity of Dow-Coming 705 oil (D.C. 705) was calculated based on 
measured transmission data. These results are shown in Fig. 7. However, the evaporation rate 
for 705 oil is higher than desired. Several new organosiloxane oils [10] have been developed 
that have much lower evaporation rates. The oil with the lowest evaporation rate is designated 
Mej in Ref. 10. To evaluate the effect of the evaporation rate on the mass of a sheet radiator 
system, consider the lifetime, t^,, of the sheet mass, M*. 


N, — P^AraD *®v*Arad 


( 6 ) 
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TABLE 1.— ABSORPTION 
COEFFICIENT FOR DOW— 


CORNING OIL 


Wavelength 

range, 

pm 

Absorption 

coefficient, 

a, 

pm 1 

0 to 2.5 

0 

2.5 to 7.0 

8.0*10“ 

7.0 to 7.7 

6.4* 10' 3 

7.7 to 8.0 

1.75* 10’ 2 

8.0 to 8.7 

1.10x10 2 

8.7 to 10.5 

1.75*10~ 2 

10.5 to 11.1 

1.10* HE 2 

11.1 to 15.4 

1.75* 10~ 2 

15.4 to 20.0 

6.4* 10" 3 

20.0 to 23.5 

5.7* 10“ 3 

23.5 to 40 

1.07* 10" 2 

40 to 100 

l.25*10' 3 

100 to 00 

0 



Figure 7.— Calculated emissivity of Dow Coming silicone oils for temperatures of 300 
and 400 K. 
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Therefore, 



( 7 ) 


where m v is the evaporation rate in kg/m 2 sec and A,^ is the total radiating area. If t L 
is the desired lifetime of the sheet radiator then the required extra fluid mass, M ST , that must 
be available to make up for evaporation loss is the following. 


M 


ST 


( \ 



V ' J 


( 8 ) 


Using the experimental results for D.C. 705 and [9,10], t* was calculated as a function 
of temperature, T. For D.C. 705 the mass loss rate, m v , was calculated using the experimental 
vapor pressure which was curve fit with the following expression [9]. 

log )0 p v = 12.31 - T = K, p v = Torr (9) 


Assuming the evaporation rate is given by Knudsen effusion, then 



( 10 ) 


where m is the mass of an oil molecule (atomic weight, of D.C. 705 = 546). Equations (4) 
and (10) were used to calculate t, for D.C. 705. However, experimental results [10] for m v 
that are curve fit with the following expression were used for Me^. 


lo g>v = 13 82 “ 


10197 

~T~ 


T = K, m = Kg/m 2 sec 


( 11 ) 


Also, experimental results [10] for p are curve fit with the following result 


p = 1221 - 0.8878 T Kg/m 3 


( 12 ) 


and are used to calculate ^ for Me^ It should be noted that m v calculated using Eq. (10) 
for Dow-Coming 704 oil in Ref. 10 is an order of magnitude greater than the experimental 
values for m v . If this same result holds true for D.C. 705 then the calculated values of m v 
using Eq. (10) may be an overestimate. As a result, the calculated lifetimes, t^ for D.C. 705 
may be greatly underestimated. 

Results of the ty calculations are shown in Fig. 8 for r = 100, 200, and 300 pm. For Mej 
the evaporative loss is so low that ^ > 10 yr for T < 440 K. Thus, the required extra fluid, 
M sx < M s for t L < 10 yr for T * 440K. For T < 400 K, M sx < 0.01 M, for t L = 10 yr. In this 
case (T < 400 K) negligible extra fluid mass is required. For D.C. 705 oil the evaporation rate 
is much higher. As a result, for t L = 10 yr the temperature, T < 300 in order for M, t < 0.1M S . 
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Figure 8. — Sheet lifetimes for Dow Coming silicone oils Me 2 and 705. 


From these results we would conclude that Mej will make a better working fluid than 
D.C. 705. However, the emissivity must also be considered. 

To calculate the emissivity of Mej the experimental spectral transmission data [10] for 
100 and 200 pm films was used. This calculation is similar to the one for D.C. 705 [7]. The 
experimental spectral transmission data is divided into wavelength bands of constant 
transmittance. From the experimental transmittance the absorption coefficient for each of the 
wavelength bands is calculated using the spectral transmittance relation for an infinite sheet 
[6,7,12]. 


= 2E 3 (a x r) 


(13) 


where a x is the absorption coefficient and ^(x) is the exponential integral. 



du 


Knowing a x the emissivity, e, can then be calculated [12]. 



e x e xb- 


dX 


(14) 


(15) 


where <j SB is the Stefan-Boltzmann constant (5.67 *10 8 W/m 2 K 4 ) and e Xb is the black body 
hemispherical spectral emissive power [12] and the spectral emittance is the following [12]. 

£ X = 1 - T x = 1 - 2E 3 (a x T) (16) 
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For the temperature range of interest (300 to 400 K), e^ b — 1 ► 0 for X < 4 pm and X > 70 pm. 
Therefore, to calculate e, it is necessary to know a A in the spectral region 4 < X < 70 pm. 
The experimental transmission data of Ref. 10 was split into 13 wavelength bands and the 
absorption coefficient for each of these bands was determined as described above. These 
results are shown in Table 1. Using this data the integral in Eq. (11) is split into 13 
wavelength bands with e A constant within each band. The emissivity can then be calculated 
using the tabulated values of 


\T'V 


= _i_f 

a T 4 A 

U SB 1 


e»<»X 


( 17 ) 


from Ref. 12. The results of this calculation for are compared to the e results for D.C. 
705 in Fig. 7. As can be seen, D.C. 705 has a higher emissivity than Me^ However, even for 
Mej, £ > 0.75 for x > 200 pm. Thus, the high absorption coefficient in the wavelength region 
8 < X < 15 pm (Table 1) results in acceptable emissivity even though the silicone oils are 
transparent in the visible wavelength region. For temperatures close to 300 K where 
evaporation loss for 705 is small, the larger values of £ make D.C. 705 more attractive as a 
working fluid than Me^ However, for T « 400 K the evaporation loss for D.C. 705 is too 
large and Mej is the more desirable working fluid. 

RESEARCH IN PROGRESS AND RESEARCH TO BE DONE 

Currently research is concentrated in two areas. First, the stability limit of the sheet flow is 
being experimentally determined. As discussed earlier, there is a limit on W/x to insure that 
holes will not form and destroy the sheet. A stability analysis is also being earned out to 
compare sheet stability predictions with the experimentally measured W/x limit. The main 
component of the experimental apparatus is a variable slit with limits of 10 pm < x < 100 pm. 

The second area being experimentally investigated is the sheet emissivity. Using the facility 
described in Ref. 7 the emissivity of D.C. 705 oil will be measured, pending program funding. 
The emissivity will be determined by measuring the enthalpy, i.e., temperature change of the 
fluid between the slit and the coalescence point. 

Collection of the sheet flow at the coalescence point is the remaining area for research. 
However, a large amount of collector work [13-15] for the liquid droplet radiator (LDR) may 
be applicable to the LSR. It may be possible to recover some of the kinetic energy of the 
sheet flow with a collector that is designed like a subsonic diffuser of a wind tunnel. The 
collector is the remaining undefined component for an LSR system. Before a prototype LSR 
can be built, the collector problem must be solved. 


CONCLUSION 

The fluid mechanics of liquid sheet flows has been extensively researched. Scaling of the 
flow has been found to satisfy the theoretical result L/W = [We/8] 1/2 . Stability of the flow 
is currently being investigated. The limit on W/x is being experimentally determined. 
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Emissivity, e, and sheet lifetime, resulting from evaporation loss have been calculated for 
two possible LSR working fluids. In the case of D.C. 705 silicone oil, its larger emissivity and 
acceptable lifetime at T * 300 K make it attractive for use at the lower temperatures. 
However, at T » 400 K the lifetime for D.C. 705 is short (< 1 yr), so that Dow-Coming Me^ 
which has C > 10 3 yr at 400 K, becomes the more desirable fluid. 

Currently an experimental measurement of the sheet emissivity is being planned. The 
remaining research problem that has not been studied is the collection of the sheet flow. 


NOMENCLATURE 

A,. Cross-sectional area of edge cylinder, m 2 

Arad Total radiating area of sheet = WL, m 2 

a x Spectral absorption coefficient, m' 1 

Ejfx) Exponential integral defined by Eq. (14) 

e b Black body hemispherical spectral emissive power, w/m 2 pm 

k Boltzmann constant, 1.38* 10 23 j/k 

L Length of sheet in flow direction, m 

M s Mass of sheet fluid, kg 

M st Required extra fluid mass to make up for evaporation loss, kg 

m Mass of oil molecule of sheet fluid, kg 

m. Evaporation rate for sheet fluid, kg/m 2 sec 

p v Vapor pressure of sheet fluid, n/m 2 

R,. Radius of curvature, m 1 

r Width of edge cylinder in x direction, m 

s Height of edge cylinder in y direction, m 

T Temperature of sheet fluid, K 

T x Spectral transmission coefficient of sheet fluid 

L Lifetime of sheet fluid, yr 

t L Desired lifetime of sheet radiator system, yr 

u e Velocity of sheet edge in x direction, m/sec 

W Width of sheet, m 2 

T W 

We Weber number for sheet flow = p — 1 
w„ Sheet flow velocity in z direction 0 

x Coordinate in plane of sheet and perpendicular to flow direction (z) 

y Coordinate perpendicular to plane of sheet and perpendicular to flow direction (z) 

z Coordinate in plane of sheet and in flow direction 
e Sheet emissivity 

e x Spectral emittance 

q Dimensionless sheet edge cylinder height = 2s/x 
A Wavelength, m 

p Density of sheet fluid, kg/m 3 

a Surface tension of sheet fluid, n/m 

a SB Stefan-Boltzmann constant, 5.67* 10 8 , w/m 2 K 4 
t Thickness of sheet, m 
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lifetime as determined by evaporation losses were made for two silicon based oils; Dow Corning 705 and Mc 2 - 
Emissivities greater than 0.75 are calculated for t > 200 pm for both oils. Lifetimes for Mc 2 arc much longer than 
lifetimes for 705. Therefore, Me 2 is the more attractive working fluid for higher temperatures (T > 400 K). 
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